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We build a realistic model of curvaton cosmology, in which the energy content is described by 
radiation, WIMP dark matter and a curvaton component. We calculate the curvature and isocur- 
vature perturbations, allowing for arbitrary initial density perturbations in all fluids, following all 
species and their perturbations from the onset of dark matter freeze-out onto well after curvaton 
decay. We provide detailed numerical evaluations as well as analytical formulae which agree well 
with the latter. We find that substantial isocurvature perturbations, as measured relatively to the 
total curvature perturbation, can be produced even if the curvaton energy density is well underdom- 
inant when it decays; high precision measurements of cosmic microwave background anisotropies 
may thus open a window on underdominant decoupled species in the pre-nucleosynthesis early Uni- 
verse. We also find that in a large part of parameter space, curvaton decay produces enough dark 
matter particles to restore WIMP annihilations, leading to the partial erasure of any pre-existing 
dark matter - radiation isocurvature perturbation. 

PACS numbers: 98.80.Cq, 98.70.Vc 



■ I. INTRODUCTION 
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■^J- , According to standard lore, the cosmological density perturbations at the origin of large scale structure were seeded 
CTS ■ at some very early time by the decay of the field that is also responsible for inflation. In the past few years, there 
has been growing interest in a variant of this cosmological scenario, in which another field, denoted the "curvaton" , 
communicates its own perturbations to matter and radiation. In the simplest version, the curvaton plays a negligible 
role in determining the dynamics of the Universe at the time of inflation, but comes to dominate the energy density in 
the radiation era, and as such, it becomes the seed of energy density perturbations through its decay. From a technical 
point of view, the curvaton scenario transforms an initially isocurvature perturbation into an adiabatic (or curvature) 
perturbation, as demonstrated in another context in Ref. [l[ and more recently in Ref. [|J (see also Ref. [H). 

In the "simplest" curvaton model, radiation is composed of a single fluid deprived of an intrinsic density perturbation 
before curvaton decay; the final perturbation is then purely adiabatic. However, as already noted in [2|], radiation 
is generically a mixture of different interacting fluids, so that if the curvaton decays after the decoupling of one of 
the constituents, an isocurvature perturbation will appear between this latter and the rest of radiation. Detailed 
^ , calculations have been provided in Refs. [fj, |7j for isocurvature perturbations generated between photons and baryons, 
between photons and neutrinos as well as between photons and dark matter and their impact on the CMB anisotropies. 
Non-gaussian signatures are also to be expected [H, Q . The implementation of a curvaton scenario in the framework 
of high energy physics and supersymmetric theories has been proposed in several studies, see for instance [Icl |. 

Since any pre-existing isocurvature perturbation between two interacting fluids is erased on a short timescale if 
the fluids share thermal equilibrium [ill] . an isocurvature perturbation can only be generated when the curvaton 
decays after the decoupling of one of these components. In this respect, the case of WIMP dark matter - photon 
isocurvature perturbation appears particularly interesting. Dark matter is indeed expected to decouple very early 
on from radiation, hence this leaves significant room for the curvaton to decay after decoupling and before big-bang 
nucleosynthesis (as required by cosmology) . This may be put in contrast with baryons and neutrinos which decouple 
just before big-bang nucleosynthesis. 

To our knowledge, there is no exhaustive survey of the curvature and isocurvature perturbations generated in 
a "realistic" curvaton model which takes into account the existence of radiation and dark matter, the influence of 
the freeze-out of this latter fluid, the influence of initial density perturbations in radiation and dark matter and the 
different decay widths of curvaton into radiation and dark matter. For instance, Ref. @ has computed the isocurvature 
and curvature perturbations under specific assumptions on the origin of dark matter, neglecting the effect of freeze- 
out and of different branching ratios. The generation of a net isocurvature fluctuation during dark matter freeze-out 
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consecutive to the presence of the curvaton has been computed in Ref. but without discussing the evolution through 
curvaton decay. In an initial density perturbation in the radiation fluid has been taken into account but the 
presence of dark matter was neglected. Ref. [13| . building on the general formalism developed in [3], has considered 
a three-fluid model of curvaton decay, albeit assuming that dark matter fully originates in the decay of the curvaton. 
Finally, Ref. 15] has presented a detailed numerical evaluation of the curvature and isocurvature perturbations in a 
three-fluid model of curvaton decay and their impact on CMB anisotropics , without however specifying the origin of 
dark matter. 

Our present objective is to make progress along these lines and to present a comprehensive model of curvaton 
- radiation and WIMP dark matter cosmology, building on the above earlier studies. We will present analytical 
and numerical calculations of the perturbations produced, following in detail the freeze-out of equilibrium of dark 
matter annihilations and curvaton decay, with non- vanishing initial perturbations in the radiation - dark matter sector 
and in the curvaton field. We focus on WIMP (neutralino) dark matter, since it offers at present one of the best 
motivated candidates for dark matter. It further allows one to set a specific framework for the thermal history of the 
different fluids, without having to rely on the details of the underlying particle physics models. The scenario, and 
its cosmological consequences, would certainly be different for axion like dark matter, whose scrutiny is postponed to 
future study. 

The paper is arranged as follows. In the following section, i.e. Sec. UU we describe the three-fluid model and 
establish the equations of motion both at the background level, see Sec. Ill Al and at the perturbed level, see Sec. IIIBI 
In Sec. IIII1 we solve these equations by approximate analytical methods. In Sec. IIV1 in order to assess the validity of 
our approximations, we investigate the evolution of the background and of the perturbations by means of numerical 
calculations. The physical interpretation of our results is also discussed at length in this section. Finally, in Sec. [V] 
we present our conclusions. Some specific technical details are also presented in three appendices at the end of the 
paper. In App. [XI we briefly recall how the gauge-invariant quantities relevant to the present work are defined and 
how their equations of motion can be derived. In App. [Bl we explain how to quantify the influence of the curvaton 
perturbations during freeze-out [our Eq. (|22"|) ] using junction conditions. In App.0 we provide an extended discussion 
of the theorem of Weinberg [Tl| on the approach to adiabaticity for the perturbations of interacting fluids in terms of 
gauge- invariant quantities. 



II. DETAILED THREE-FLUID MODEL 



We are interested in computing the curvature and isocurvature perturbations produced in a three-fluid model 
composed of radiation (denoted by the subscript 7 in what follows), a curvaton component (denoted a) and dark 
matter (denoted by %). The dark matter particles share thermal equilibrium with radiation until their freeze-out 
of equilibrium after having turned non-relativistic. We denote by Xf the ratio of mass m x to temperature at which 
freeze-out occurs: Xf = m x /Tf. In the particular case of neutralino dark matter, one generically finds Xf ~ 20. Unless 
otherwise noted, we use Xf = 20 as a fiducial value in what follows. 

We further assume that the curvaton is very weakly coupled to the visible sector (radiation + dark matter) and that 
it decays into this sector after dark matter freeze-out. Were it to decay earlier, there would be no final isocurvature 
perturbation between x an d 7, in virtue of the theorem of Weinberg on the approach to adiabaticity of interacting 
fluids [11]. This point is addressed with greater scrutiny in Appendix [UJ where it is shown, in particular, to be in 
excellent agreement with numerical computations. If the curvaton were to couple to the visible sector so that it had 
at some point shared thermal equilibrium, any pre-existing isocurvature perturbation between itself and radiation or 
dark matter would likewise have been erased. Hence the above assumption is one of non-triviality. It is furthermore 
consistent, in the sense that a late decay is a generic consequence of a very weak coupling; moduli fields, in particular, 
emerge as natural candidates in this curvaton scenario. 

We denote by T a the total decay width of the curvaton, r cr7 and T ax its partial decay widths into radiation and 
dark matter, respectively. Decay of the curvaton occurs when T a ss H , hence our previous assumption corresponds 
to To- <C r xx | f , where the latter quantity r xx | f = n x \ f (a xx v) denotes the x ~ X annihilation rate at freeze-out. The 
annihilation cross-section is related to the time of freeze-out by the standard formula [16| : 

with mpi = G -1 / 2 ~ 1.2 x 10 19 GeV, g — 2 and g* ~ 230. In the following, we describe how the previous physical 
situations can be modeled at the background and perturbative level. 
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A. Equations for background quantities 

In the non-relativistic era, and in the absence of dark matter production from curvaton decay, the number density 
of dark matter particles decreases due to expansion and annihilation, see e.g. [16| : 



An 



x 



+ 3Hn x = - {<7 xx v) (n 2 x - n\ >eq ) , (2) 



where n x . cq is the density at thermal equilibrium defined by 

-oo ( E 2 _ m 2\ 1/2 E n „ o,, 

where <? is the number of spin states. The last equality assumes that the dark matter particles have become non- 
relativistic, i.e. x = m x /T S> 1. Since p x ~ m x n x in the non-relativistic era, the equation of evolution for p x takes a 
similar form. In order to take into account a source term arising from curvaton decay, it suffices to add a term Y ax p a 
in the r.h.s. of the above equation. One thus arrives at 

^ + zhpx = (pI - p 2 x , q ) + r.x/v ■ (4) 

Dark matter annihilations produce radiation, consequently the evolution of the energy density in the radiation fluid 
is described by a similar equation: 

^ + U*(h = +^ (Pi ~ P'U + ■ (5) 

Finally, the energy density of the curvaton component decreases through expansion and decay: 

+ 3Hp a = - T a p a , (6) 

where, since r CT = T ax + r CT7 , one has conservation of the total energy density. Note that we treat the curvaton as a 
non-relativistic fluid or as a coherent scalar field oscillating in a quadratic potential. 

These equations, as well as the Friedmann equations, can be recast in dimensionless forms in terms of the density 
parameters f2 7 , fl x , ft a , the Hubble factor H, and the e— fold number N (or equivalently the parameter x = m x /T), 
which is defined in terms of the scale factor a by diV = din a = dlnir: 



dA H + x 7 Stt m [ * ' [n 



x 

1>^-<M<»- - D - -^ — (ii; <>;.,„) , (8) 



dn -r _ a r »un ro n i 3 ( CT xx u ) m Pi g 



We assume that radiation is always in a state of thermal equilibrium, so that the parameter x which enters the formula 
for jO x . oq can be written in terms of the radiation energy density, hence H and £! 7 : 

, x = /TrV J7T \ 1/4 _ (n) 

For the sake of simplicity, we ignore any temperature dependence of the function g„ . 



B. Equations for perturbed quantities 



We now turn to the description of the equations of motion for the perturbed quantities. In Appendix \^ it is 
shown how to find the corresponding expressions from the general perturbation of a fully covariant expression for 
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the annihilation and decay terms. Using these results, one can write the evolution equations for the gauge invariant 
dimensionless quantities A( Q ) = Ap( Q )/p( Q ) with a = 7,x> f, and <&: 

x [-* (n> - f^ oq ) - 2 (O^A x - nl cq A cq ) + (Q 2 X nl cq ) A x ] , (12) 



i (0 CT A ff + 7 A 7 + O x A x ) - <E ( 3 + ^ t M , (14) 



dTV 2 v " " ' 1 x V // 

d$ 1 

= -$ 

AN 2 



-$--(n CT A ff + n 7 A 7 + n x A x ) , (15) 



where we have introduced the short-hand notation 

1/3 m\ . 

AoqS 4 [2 + t) A ^ (16) 

This formula follows from Eq. ([3]), writing the temperature fluctuations in terms of Ap 7 according to Eq. (fTTj) . Note 
also that, throughout this study, we assume that (<r xx v) is a constant, independent of temperature, meaning that 
annihilations only occur through s— wave interactions. 

It is also interesting to re-write these equations in the following form: 



diV V " ' 3~V H n Y { x a> H n Y 8ttto v fi 



X JJ "X o/i//t x 

$ (fi* - f£ „) - 2 (fi x A x - fi x , eq A eq ) + (Ql - fi» ) A x ] , (17) 



X X<cq/ V X X X,eq e( U ~ \ X X,cq 

W + r H Cl-y a,+ H 87rm Y fl 



p-y _t J. J L-y U / l / < 6^ d icy 

x [* ("? - n U) + 2 K a > - "L>AJ - ("5 - A J . ( ls ) 
3 ^(-* + H =-*H^- do) 

When the r.h.s. of the previous expressions vanish and/or are negligible, the equation of state of each component does 
not evolve in time, and one finds that the quantities (where a dot means a derivative with respect to cosmic time) 

« a) = -*-H^ + , (20) 

P{a) 3|_1+w (q) J 

are conserved, a well-known result for isolated fluids. In this expression, = P( a )/P(a) is the equation of state 
parameter for fluid a. 

III. CURVATURE AND ISOCURVATURE PERTURBATIONS: ANALYTICAL CALCULATION 

We now compute the isocurvature and curvature perturbations that arise in this three-fluid model using analytical 
approximations for the freeze-out of dark matter and for the decay of the curvaton. In the next section, we compare 
the results of numerical evaluations of the above equations to these analytical calculations. 

Since we assume that the curvaton never interacts with radiation nor dark matter until it decays, Co- is conserved 
(provided that the equation of state of a does not change, which we assume to be the case). Until dark-matter 
freeze-out, radiation and dark matter share thermal equilibrium, hence £ 7 = £ x as shown by Weinberg (llj |. see also 

Appendix [C] Therefore, we can set the initial conditions just prior to dark matter freeze-out, and write Q ■ , and 
C 7 J) = Cx the initial gauge-invariant Q a \ quantities in curvaton, radiation and dark matter at that time as defined in 
Eq. (12TJ)) . We use standard notations for the isocurvature modes, namely: 

S xl = 3 (( x — C7) , S a ™ = 3 (C CT — C7) • (21) 
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FIG. 1: Comparison of numerical (red solid line) and analytical (dotted blue line) calculations of the isocurvature perturbation 
S xl generated by the curvaton during the freeze-out of \ ~ X annihilations, before curvaton decay. 



Hence, initially S x \ = but is generically non-zero. 

In the presence of a curvaton, a net isocurvature perturbation S XS) may be generated at two different epochs: during 
dark matter freeze-out, and at curvaton decay. Lyth and Wands [9j have shown indeed that S X7 ^ after freeze-out 
since the very definition of freeze-out, T xx = H defines a hypersurface which does not necessarily coincide with 
the hypersurface of uniform radiation energy density. In other words, on this latter hypersurface, dark matter and 
radiation do not present density perturbations prior to freeze-out, yet freeze-out does not occur at the same time at 
all space points as a result of the non-uniformity of the Hubble factor on this hypersurface. In Appendix [B] we derive 
the resulting isocurvature perturbation using the properties of junction conditions in general relativity. We obtain, in 
perfect agreement with Lyth and Wands [§] : 

c>! _ c<t , (<*f - 3 ) f (r<l _ t <<\ _ c<! , 1 (of - 3) n>* <f 

^ ^ 2 (a f - 2) + fi> f [ ^ ^ ' + 32(a f -2) + n> f CT7 ' [ ' 

where quantities denoted by a superscript >f (resp. <f) are evaluated just after (resp. just before) freeze-out. Notice 
that we could also replace C x * by £5* because of the Weinberg theorem as discussed before. The variable a is related 
to the scaling of the annihilation rate with temperature; in practice, at = Xf + 3/2 ~ 20. Freeze-out of annihilations 
does not affect radiation nor the curvaton, hence one can safely assume: £> f = £< f and £> f = £< f , as long as curvaton 
decay has not yet occured. Therefore, a net isocurvature perturbation S xl is generated from the pre-existing S ai 
curvaton - radiation isocurvature perturbation: 

_ (of-3)n>* m 

b ^ ~ 2(a f -2) + n> f6 ^ ■ [ ' 

As expected, S x ^ vanishes in the limit fl^ 1 — > 0. In Fig. [TJ we compare the above analytical approximation with 
detailed numerical integration of the equations for the perturbed quantities that were presented in the previous section. 
There is overall agreement although the two curves tend to differ by a factor of a few at small values of f2^ f , the 
energy density parameter in the curvaton at the time of freeze-out. 

During curvaton decay, a net isocurvature perturbation 5 X7 may also be produced if either dark matter or radiation 
(but not both) inherits (part or all of) the curvaton perturbation. This isocurvature perturbation can be calculated 
as follows. In the time interval following dark matter freeze-out and preceding curvaton decay, the various Q a ) are 
separately conserved since the various fluids do not interact with each other. In order to compute the evolution of 
the curvature perturbations through curvaton decay, one may use the total £, namely 

C = E — = rCr + < 7 + (1 - r - s) ( x , (24) 

(a) 9 



6 



where 

30 CT + 3fl x + 4Q 7 30 CT + sn x + 4ft 7 ' v ; 

Although £ evolves in time due to the evolution of r and s, it is conserved through decay if one assumes curvaton 
decay to be instantaneous. Hence one can relate the quantities immediately after decay, denoted by the superscript 
>d, to those immediately before decay, denoted <d, through the relation: 

r <d (< d + s <d C< d + (1 - r <d - s <d ) (< d = t >d (> d + (1 - < >d )C> d , (26) 

with 

We may further assume that the energy density contained in dark matter immediately before and immediately after 
decay is negligible. Curvaton decay must indeed proceed before big-bang nucleosynthesis, and dark matter domination 
well after. Then we may set t >d ~ 1, (1 — r <d — s <d ) ~ 0, so that: 



C> d = (1 - r <d )C< d + r <d C< d - C 7 d + ^ <d S< 7 d . (28) 



This result for £ 7 had been obtained previously in two-fluid models of the curvaton scenario. If one removes the 
assumption t >d ~ 1 (which implicitly implies £? 7 = Y ai /Y a ~ 1), it is possible to show that, in the limit £? 7 — > 0, 
one finds (> d — > (< d independently of r <d , as it should. Moreover, the coefficient r <d oc Yl< d and, therefore, in 
the absence of a curvaton, the quantity £ 7 is unaffected as expected. Finally, if there is no pre-existing isocurvature 
perturbation between the curvaton and radiation, £ 7 is also left unchanged as the perturbations share a common 
origin. 

Let us also remark that, since (< d ~ C^"' — C^ f — Cr an d similarly (< d ~ Q 1 ', the above equation allows to express 
the final perturbation in radiation in terms of the initial perturbations. This will be done further below. 

In order to evaluate the final curvature perturbation in dark matter, it is convenient to construct a compound fluid 
that comprises the dark matter and a fraction B x = Y ax /Y a of the curvaton fluid such that the source term coming 
from curvaton decay vanishes in its equation of energy conservation (l3j : 

Pcomp = p x + B x p a . (29) 

For this compound fluid, the curvature perturbation ( comp is conserved, with: 

^ ^ Qv -— B x £l a * £l x -\- B x Yl a ^ 

The conservation of Ccomp comes from the fact that the source terms also cancel out at the perturbed level. Notice 
that this is possible only because the curvaton and the dark matter have the same equation of state; in particular, 
one could not construct a simple compound fluid made of radiation and the non-relativistic curvaton. Then, writing 

Comp = Comp and = leads to: 

yi >i r n >f i r n >! 

^X o>f _L R 0>f 0>* _1_ H > f^ cr ^X 1 0>f _L R 0>f a X ' ^ > 

\L X -f- d x \Lo \L X -f- r> x \L(j o\L x -\- £> x \Lcr 

This equation has a structure very similar to Eq. I|28j) . It indicates that dark matter perturbations inherit a contri- 
bution from curvaton perturbations only if the branching ratio B x , the curvaton contribution before decay and the 
pre-existing S^ d do not vanish, as is intuitively clear. 

Combining Eqs. (|31[) and (|28p . one may thus write the isocurvature perturbation S X7 generated through curvaton 
decay as a function of the pre-existing S x -y, S xa and/or S^y. 

R >f 

c>d _ c<d I 1J X iL a o>f „<d c<d 

a xy ~ ^x-r + n >! + Bx ^i>i ^°x ~ r a <n ■ y 6l > 



In the r.h.s. of the above equation, terms can be evaluated at any time in the interval from post-freeze-out to pre-decay, 
except r which must be calculated immediately before decay. 
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To summarize, the final radiation curvature perturbation can be expressed as 

<*) = Cf + ±r«S« . (33) 

This final Q > is in fact nothing but C7 given by Eq. ([28]) because the perturbations in the radiation fluid (and in 
the curvaton) are not affected by the freeze-out. On the other hand, using Eqs. (I31|) and (|28|) . the final dark matter 
curvature perturbations can be written as the sum of two contributions, one from the freeze-out and one from the 
curvaton decay, namely 

(f) 1 (a f -3)n>* 1 B x Sl? 

^ x ^ x 3 2(a f -2) + ft> f CT7 3f2> f +B x ft> f CT * 

1 (a f -3)n>* »> f <f 1 S x n>* <f 

" ^ x 3 2 (a f - 2) + n> { Q> ! + B x fl> { + 3 n> ! + B x n> ! yo ' 



Therefore, combining Eqs. (|33[) and (|34[) . our result for the final isocurvature perturbation in terms of the initial 
isocurvature perturbation can be written as 



(a f -3)n>< ^> f S x O> f 



2(a f - 2) + n> f ^> f + B x n> ! n> { + B x n> { 



XI 

\_*\u,i - *) T a iff » X T -u x " CT 

(i) 



S« , (35) 



where we have used the results: £v f = C7 and S 1 ^ = = S a \. Assuming radiation domination all throughout 
freeze-out, the energy densities just after freeze-out are given by [see Ref. [IB], especially Eq. (5.45)] 1 : 

J2>* = fl^ — , tt> { ~ 1.67 x 10" 3 % /2 e- Xf . (36) 

It proves interesting to discuss various limits of Eq. (f35|) . Consider for instance the case B x tt> { ^> f2> £ . Physically, 

this means that curvaton decay will produce more dark matter than existed prior to decay, hence dark matter entirely 

inherits the curvature perturbation of the curvaton: <~ ffl and consequently, S X j ~ (1 — r <d )Sa\. In the limit in 

which the curvaton density parameter is negligible at the time of decay, r <d <C 1, only a negligible fraction of radiation 

has been produced during curvaton decay, hence radiation conserves its initial curvature perturbation: ~ Q 1 ', 

({) 

and therefore the isocurvature perturbation S x <y is maximal. We find it important to stress that in this limit, a large 
isocurvature perturbation may be generated during curvaton decay, even if the curvaton contributes a negligible amount 
to the total energy density at the time of its decay. This result comes about because the energy density of dark matter 
in the very early Universe, prior to big-bang nucleosynthesis, is itself very small, hence it can be strongly affected by 
the decay of underdominant species. 

In the opposite limit, r <d ~ 1, meaning Q< d ~ 1, the decay of the curvaton regenerates both radiation and dark 
matter, hence the final isocurvature perturbation turns out negligible. 

It may also be that the branching ratio B x of curvaton to dark matter is negligibly small, or even zero, for instance 
if the decay is kinematically forbidden. In this case, a net isocurvature perturbation may arise from two different 
sources: freeze-out and/or transfer of a perturbations to radiation, as can be seen in Eq. (|35p . The former depends 
on the magnitude of the curvaton density parameter at the time of freeze-out: it is negligible if Q> 1 <C 1, but tends 
to « 50 % as — * 1. The latter increases with r <d , as it should. In particular, S x j is maximal if fijT* <c 1, so that 

Q; ~ £ x \ but r <d ~ 1, so that Q ~ Q- , i.e. all of the curvaton perturbations is transfered to radiation but none 
to dark matter. 

These limits will be encountered in the following Section, in which we compare the above analytical results to 
numerical integration of the equations derived in the previous section. 

IV. CURVATURE AND ISOCURVATURE PERTURBATIONS: NUMERICAL EVALUATION AND 

DISCUSSION 

The above analytical calculations assume instantaneous dark matter freeze-out and instantaneous curvaton decay. 
In order to assess the validity of the previous results, we have integrated numerically the equations of motion given 



1 We compute using the variable Y&0 defined in this equation and assuming radiation domination at freeze-out; this notably guarantees 

>f 
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in Section [H] starting before freeze-out at x\ = 3, until well after curvaton decay. In general, the previous analytical 
formulae are found to be accurate but they tend to overestimate the final isocurvature perturbation in a substantial 
fraction of parameter space, by a factor sometimes as large as an order of magnitude. The source of this error lies in 
the neglect in the above analytical calculations of the coupling between dark matter and radiation after freeze-out. 
Indeed, when a significant amount of dark matter particles is produced by curvaton decay, annihilations may become 
effective again, which leads to a erasure of the pre-existing S xl isocurvature perturbation. This will be detailed 
and quantified in what follows. The numerical integration also allows to follow the evolution of the curvature and 
isocurvature perturbations in the different components in parallel with their abundance, which as we will show, offers 
a sound and intuitive understanding of their behavior. 

A. Background and perturbations evolution 

The parameters that control this evolution are obviously: Via , which controls the relative magnitude of the curvaton 
energy density at decay and at freeze-out, IV (as measured relative to the fixed quantity r xx |f ) which controls the 
time of decay, the branching ratio B x = T tTX /T tT , which controls the amount of dark matter produced relative to 

radiation during curvaton decay, and C7 /C^ > the- ratio of initial density perturbations in radiation (or dark matter) 
and the curvaton. However, the value of CyVC^ on ly influences the final total curvature perturbation ffl [or Qp] 
as expressed relatively to Q 1 , see Eq. (|3"5)) . The isocurvature perturbation S x ]/Sal, i.e. measured relatively to the 
initial isocurvature perturbation, does not depend on this parameter, see Eq. (|35p . We may therefore keep CtVCt^ 
to a fixed value (here, 0.31) without affecting the generality of our results. 

In principle, one could determine these various parameters if the couplings of the curvaton to the differents sectors 
of the theory as well as the early dynamics of the curvaton and the inflaton and their perturbed components during 
the inflationary era were specified. This however lies well beyond the goals of the present study. Furthermore, our 
objective is to present a complete survey of the various alternative cosmological effects, hence we have chosen to select 
reasonable values for these parameters in order to bring to light the most relevant consequences. 

The time of freeze-out, encoded in X{, also affects the discussion since it directly determines the dark matter energy 
density. We will keep this value fixed to the generic X{ — 20 in most of our discussion, but discuss its influence, and 
show one particular case of the evolution in the case X{ — 30. We have also decided to keepthe time of radiation- 
matter equality free although it should of course be fixed to T equ ~ 5.5floh 2 eV (x eqn ~ 10 11 ) It is always possible 
to achieve the right dark matter abundance by tuning the above parameters, but this would be obtained at the price 
of unnecessary complexity. However, we set the time of curvaton decay to take place before big-bang nucleosynthesis, 
which starts at x BBN ~ 10 5 (m x /100 GeV). 

In order to survey the various effects in this extended parameter space, we plot in a first series of figures typical cases 
of the evolution of the background and perturbed variables as a function of e— foldings of the scale factor (actually 

x = m x /T) for various values of Q$ and B x . 

In Figure^ the evolution of fl x (upper panel, dashed line), f2 7 (upper panel, solid line), fl a (upper panel, dotted 

line), Ct/Co- (middle panel, solid line), Cx/Co- (middle panel, dashed line), and S xl /Sa\ (lower panel), is displayed 
as a function of m x /T for an initial value Qa = 3 x 10~ 7 , which corresponds to = 2.0 x 10~ 6 [sec Eq. (|36p ] 
and to a "small" value of Sl< d = 2.7 x 10~ 4 , that is to say a small value of r <d . The branching ratio is taken to be 
B x — 10~ 3 ; the curvaton decay width is set to T x — 10 _6 r xx |f , which corresponds to a time of decay Xd = 4.7 x 10 3 
(defined as the time at which H = T x ), or, in terms of temperature, T >d ~ 21 MeV (assuming m x = 100 GeV). 

As shown in the upper panel of this figure, the Universe remains radiation dominated at all times before radiation- 
dark matter equality. The curvaton energy density increases relative to that of radiation in proportion to a(t) oc x, 
which is characteristic of non-relativistic matter, then vanishes exponentially as the curvaton decays. The dark 
matter energy density parameter first decreases exponentially due to annihilations in the non-relativistic regime, then 
stagnates when freeze-out occurs at Xf = 20 and finally increases linearly with a, just as the curvaton, in the post- 
freeze-out era. The curvaton decay produces a significant number of dark matter particles at x ~ 0.5 — 1.0 x 10 4 , 
which translates in a faster increase of Q x in this region. 

The middle panel in this figure illustrates the concomittant evolution of the curvature perturbations Q a in these 
components. The spikes that appear for CilQx an d Cxi Or m the pre-freeze-out era come about because the curvature 
perturbations involve the inverses of the time derivatives of p 7 and p x , see Eq. ([20]) . which may become singular locally 
when production/annihilation terms are present in the equations of motion ; this, nevertheless, does not signify that 
perturbation becomes inapplicable but rather that the spatial hypersurfaces are ill-defined at that point [14j . 

Since the curvaton and dark matter are so much underabundant with respect to radiation, this latter fluid is 
isolated in a good approximation, hence its curvature perturbation £ 7 remains conserved, as seen in the middle panel. 
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FIG. 2: Evolution of background quantities (upper panel) fl x (dashed line), fi 7 (solid line), f2 CT (dotted line), curvature 
perturbations (middle panel) Cy/C^ (solid line) and Cx/C^ (dashed line), and the isocurvature perturbation S^/S^y (lower 
panel). Parameters are: B x = 10" 3 , IV = 10 _c T xx | f , and Q$ — 3 x 10" 7 corresponding to Qa d = 2.7 x 10~ 4 . 




FIG. 3: Same as Fig.[2]for the following parameters: B x = 1(T 6 , r CT = lCT 6 r xx | f and = 1CT 2 corresponding to f^ d = 0.65. 



In contrast, the dark matter perturbation deviates from its initial value at x ~ 10 3 — 10 4 as it acquires the curvaton 
perturbations through the decay of this latter. The radiation-dark matter isocurvature perturbation then becomes 
non- vanishing. 

The physical situation depicted in this figure is particularly interesting because it confirms that a significant isocur- 
vature perturbation S X j /S&\ ~ 90 % is generated even though the curvaton field never contributes significantly to the 
total energy density content. This stems from the very small energy density of the dark matter component at the 
time of curvaton decay relative to the amount of energy transfer from curvaton to dark matter at that time (B x Q< d ). 
Equivalently, in terms of these quantities evaluated at the time of freeze-out: B x £l> { ~2x 10~ 9 > ~3x 10~ 10 , 
see Eq. (|36p. Accordingly, one can see ( x deviate from its initial value at the time of curvaton decay in Fig. [2] and the 
value of S xl evolves simultaneously. 
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FIG. 4: Same as Fig. [2] for the following parameters: B x = 10 7 , r CT = 10 6 F xx |f and C$ = 3 x 10 1 corresponding to 
flv d = 0.66. 

Figure [3] shows a similar case, but with a higher initial curvaton energy density, flS = 10~ 2 , and a smaller branching 
ratio to dark matter (B x = 10~ 6 ). This implies Q> ! = 6.7 x 10~ 2 and Q< d = 0.65. 

The evolution is similar to that seen in the previous figure, although in the present case, the curvaton comes to 
dominate the energy density before it decays, which implies that fi 7 decreases somewhat during this period, before 
being regenerated during curvaton decay. The middle panel shows that the curvature perturbation of dark matter 
tends toward that of the curvaton well before decay at x ~ 10 2 — 10 3 . Indeed, the time at which £ x departs from its 
initial value is set by matching the amount of dark matter produced by curvaton decay with that initially present. 
Furthermore the energy density in the curvaton is larger here than in the previous figure, hence this equality occurs 
earlier. Since the limit B x fl> 1 > f2> f remains valid, a net isocurvature sets up, but gets erased as radiation is in turn 
regenerated by curvaton decay (r <d ~ 1 in this case). 

In Fig. [U Cla = 3 x 10 _1 so that the curvaton dominates the energy density from x ~ 10 to x ~ 6 x 10 4 , at which 
time it decays. Correspondingly, the dark matter energy density parameter remains constant in the post-freeze era, 
until it becomes sourced by curvaton decay at x ~ 10 3 ; the radiation energy density parameter decreases in proportion 
to l/a(t) ocl/x until curvaton decay. As increases, its influence during freeze-out becomes noticeable, see Fig. [4] 

where Qa =3x 10 _1 and B x — 10~ 7 . The isocurvature perturbation S xl is then generated at freeze-out due to the 
effect discussed by Lyth & Wands @. Closer to curvaton decay, the dark matter perturbation tends further toward 
that of the curvaton. However, as in the previous example, B x il> f > fi^ f , r <d = fl< d ~ 1 so that the isocurvature 
perturbation gets erased as the radiation perturbation catches up with that of dark matter. 

Finally, Fig.Oprovides an example of the case in which the decay of curvaton to dark matter is forbidden, i.e. B x = 0. 
Since Q^ d ~ 1, all of radiation is produced during the decay. However, dark matter nevertheless inherits part of the 
curvaton perturbation through the influence of this latter on the freeze-out, as is apparent at times x ~ 30 — 10 3 
in Fig. [5] Therefore the final S xl isocurvature perturbation is not maximal, but of the order of 70%, as measured 
relatively to the initial S ai . 

B. Final isocurvature perturbation 

In a second series of figures, we explore the dependence of S x \/Sa\ on the initial abundance of the curvaton (or, 
equivalently, its abundance at the time of decay, Q^ d ) for different values of the branching ratio B x of curvaton to 
dark matter, and compare the numerical results to the analytical formulae. 

In Fig. [6l the branching ratio is chosen as B x = 0.1. Let us describe this plot in more details. For very small values 

of fl^ , there is no isocurvature perturbation (as expected). This is also the case for large values of Qa , i-e. when 
dark matter and radiation mostly originate in curvaton decay and, hence, have a common origin. Moreover, from our 
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FIG. 6: Final isocurvature perturbation S xl expressed relatively to the initial S a -f isocurvature perturbation (left y— axis), as 
a function of numerical calculation (solid line), the analytical approximation given by Eq. (|35p (dashed line) and a refined 
approximation taking into account the late time annihilation due to the curvaton decay (dotted-dashed line). The dotted line 
indicates the corresponding value of (right y— axis) for convenience. The parameters are chosen to be: F CT /r xx | t = 10 -6 
and B x = 1CT 1 . 



previous discussion of the analytical calculations, one expects that S X j will be maximal as long as B x fl> { > 17> f but 
(1 — B x )fl< d <C fl^ d ~ 1; this region corresponds to the peak plateau of the dashed line in Fig. [6] Indeed, combining 
Eqs. ([55]). the condition B x fl> 1 > fJ> f reads: 

B x n^ > 5 x 1Q- Z x] ,2 e- Xi , (37) 

with xV" 1 exp(— X{ ) ~ 0.92 x 10~ 8 for X{ = 20 while the condition (1 — B x )fl< d <C 1 can be written as: 

(l-IJ x )nW V2 oxr 1 , (38) 

V 1 xxlf / 
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since one assumes radiation domination all throughout and x\ — 3. Recall that the ratio r CT /r xx |f represents the ratio 
of the Hubble factors at curvaton decay and at freeze-out. 

As an example, consider the case depicted in Fig. [6] Xf = 20, r CT /r xx | f = 10~ 6 and B x = 10 _1 . Condition (|37|) 

shows that significant isocurvature perturbation can be produced when fi^ > 4.6 x 1CP 10 , which agrees well with the 

results shown in Fig. [6] Condition (|38|) further limits this region to fl^ < 1.7 x 10" 4 , which, again, agrees with the 
limit of the dashed line (analytical calculation). 

Despite the above considerations, one notices in Fig. [6] that, when S X j/Sa^ is significant and not too small, 
there is a mismatch between the analytical and numerical results by a factor of a few to an order of magnitude over 
most of the above parameter space. The reason of this discrepancy has been alluded to earlier; it lies in the neglect of 
the annihilation term in the post-freeze-out era. If a sufficient number of dark matter particles are produced during 
curvaton decay, this term may become significant, even at values x S> Xf, and regenerates some form of coupling 
between dark matter and radiation, leading to the partial or nearly complete erasure of the isocurvature perturbation. 
One can actually observe a typical example of the evolution of S xl under the influence of regenerated dark matter 
annihilations in Fig. [5] at times x ~ 10 36 — > 10 4 : the isocurvature perturbation decreases slightly to achieve its final 
value during this time interval; the impact of annihilations on Sl x is not apparent in this figure but becomes clear 
if one plots Vt x /a cx a 3 n x . Note that dark matter annihilations at around the time of big-bang nucleosynthesis may 
have interesting phenomenological consequences and distinct astrophysical signatures, see notably Ref. [17| . 

One can assess the range of parameters in which the annihilation term becomes significant, hence the analytical 
calculation inaccurate, as follows. If curvaton decay is instantaneous and all of dark matter originates in this decay 
(so as to produce a substantial 5' X7 ), one has: f2 x d ~ B x fl< d . After curvaton decay, one must compare the influence 
of the terms scaling with {<J xx v) in Eq. (|12p with the terms due to expansion. Since A x and $ are of the same 
order of magnitude, it suffices to compare the prefactors; furthermore, ^ x ,eq ^ when x 3> Xf, so that, finally, the 
analytical calculation should be valid provided: 

\ xx 1 P1 ^v < 1 . 39 

8ttto x x y ' 

Of course, it is understood that all terms in this equation should be evaluated immediately after curvaton decay. 
Since Q x — n x m x /p, using the Friedmann equation one can rewrite the above constraint in the more intuitive way: 

T = < 1 , (40) 

which indeed expresses the fact that dark matter annihilations do not occur since the annihilation rate remains smaller 
than the Hubble rate. Note also that, in this sudden decay approximation, annihilations cannot occur in the interval 
separating freeze-out from curvaton decay. 

Alternatively, one can discuss the case in which curvaton decay is not instantaneous, assuming here as well that 
f2 x d ~ B x il< d . One must now compare the influence of the terms scaling with (c xx v) in Eq. (fT2"|) with the term due 
to curvaton decay, which leads to the following condition: 

3ff(<7 xx u)m^ 2 T 

— — n x «B x -n a , (4i) 

at all times before curvaton decay. By considering the time evolution of the various terms, one can convince oneself 
that the l.h.s. increases faster with x than the r.h.s. Let us discuss this last point in greater detail. The r.h.s. scales 
as cx a; 3 since fl a cx x and H cx x~ 2 in a radiation dominated phase. After freeze-out, but before curvaton decay, one 
might think that, in Eq. ([7]), only the term 17 x f2 7 matters. In this case, f2 x cx x and the l.h.s. is constant, i.e. does 
not increase faster than the r.h.s. However, the influence of curvaton decay is felt before the time H — T a so that the 
term proportional to £l a in Eq. cannot always be ignored. When this term dominates, fi x cx x A hence the l.h.s. 
indeed increases faster than the r.h.s. This scaling will last as long as the annihilation term remains negligible. When 
this latter becomes important, the evolution becomes such that the annihilation term and the term proportional to 
tt a balance each other, thus scaling in the same way, which leads to fl x cx x 5 ^ 2 after the period where f2 x cx x 4 and 
until H = T a 2 . These considerations justify the fact that the l.h.s. increases faster than the r.h.s.; hence it suffices to 



In fact, neglecting the term proportional to X;e q, Eq. 10 can be solved exactly. The solution reads 



)>f -au 2 h(v) + aw 5 / 2 T 1 / 2 / (f) - /3« 2 i^i(t>) - /3u 5 / 2 T 1 / 2 Ko(«) 
[ah(u) + pKxiu)] 
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FIG. 7: Same as Fig. H but for B x = 0. 



evaluate the inequality (|41[) at the time of curvaton decay: if it is then satisfied, it will have always been verified prior 
to decay. After decay, of course, curvaton decay does not source anymore the evolution of f2 x and all perturbations 
stay constant, as seen in previous figures. Since the time of curvaton decay is set by H ~ r CT , and fl> d ~ i? x f2< d , 
the previous condition of validity of the analytical calculation can be simplified down to the same condition (|4"0")) than 
obtained in the sudden decay approximation. In terms of our parameters, this reads: 

/ r \i/2 

B x n$ —2- « 5.0 x lO-^V"* . (43) 
V i xx|f / 

Hence, the r.h.s. becomes 1.0 x 10~ n for Xf — 20. For the example considered above, namely Y a /T xx ^ = 10~ 6 and 

B x = 10 _1 , condition (|43|) shows that the isocurvature perturbation remains maximal as long as: ^ 4.6 x 10~ 7 , 
in good agreement with Fig. [6] 

In order to model the effect of regenerated annihilations on the perturbations, one can proceed as follows. One can 
study the evolution of £ x and £ 7 after curvaton decay, assumed to be instantaneous for simplicity. In this case, one 
can safely neglect the first term in Eq. ([7]) since it is proportional to £l a and the last term proportional to ^ x ,°q as 
before. If one further considers that radiation dominates (hence Q 7 ~ 1), then Eq. ([7]) can integrated exactly and the 
solution reads 

o - o>d (x/x d ) 

?[(x/x d )-l] + (x/x A )' (44j 

where x d is the time of curvaton decay. Then, making the same assumptions as before in Eq. (|17p and considering 
that $ remains constant (which is well- verified numerically), one obtains 

A x + $ = [A>* + ^} ix/Xd ^ . (45) 

T [(x/xd) - 1] + (X/Xd) 

Using this solution and the fact that A 7 ~ —2$ is constant, one can estimate the amplitude of the isocurvature 
perturbations well after the curvaton decay. Straightforward manipulations leads to 

5 (f) = S >d H $ >d (46) 

°xi i + T X7 2(T + 1) { ' 



where / and K are modified Bessel functions and u = x/xf, v = 2Y 1 / 2 « 1 / 2 [T being defined in Eq. ( 1401 )1 . The quantities a and f3 are 
two constants that should be chosen such that fl x = Q x e when u = 1 and such that the derivative of fl x has the value indicated by 
Eq. ||40[ |. for instance also at u = 1. Using the asymptotic expansion of the Bessel functions for large arguments, it is easy to prove that 
Q x oc x 5 / 2 as claimed in the text. 
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FIG. 8: Same as Fig. H but for x f = 30. 



This result indicates that the final isocurvature perturbation is suppressed by a factor (1 + T) which becomes all the 
more important as annihilations are effective at curvaton decay. This is reminiscent of the Weinberg theorem of the 
approach to adiabaticity (i.e. the erasure of isocurvature perturbations) of interacting fluids (here, dark matter and 
radiation through annihilations), although the physical conditions are here quite different, see Appendix [C] 

From our numerical experiments for different values of the parameters, we find that this reduction factor reproduces 
well the numerical results. For instance, in Fig. [5J the suppression factor manifests itself by a reduction of Sxj for 

^o- ^ 4.6 x 10 -7 , as indicated by the dotted-dashed curve which agrees well with the numerical calculations. However, 
the presence or functional form of the second term in the r.h.s. of Eq. (|4"6"|) is not confirmed by numerical experiments 
in which S^y = 0; one finds, in particular, S X y — in this case. This error is attributable to the assumption of 
instantaneous curvaton decay, as we have checked numerically. Nevertheless, even if this term is included in the 
formula (as done in the comparison to the numerical results in Figs. I6I8| ). the agreement remains quite reasonable. 

Let us now describe Fig. [7] It shows an altogether different scenario, in which the branching ratio B x = 0, i.e. the 
curvaton decays only into radiation. In this case, significant isocurvature is produced when £!^ d ~ 1, as discussed 
before. In this case, the analytical calculation always agrees quite well with the numerical results. The isocurvature 
perturbation tends to decrease as f2^ exceeds ~ 10~ 2 because the curvaton energy density at the time of freeze-out 
~ 60^ for Xf — 20] then becomes significant. As a consequence, the dark matter curvature perturbation is 
affected by that of the curvaton, hence the final dark matter - radiation isocurvature is reduced. 

Let us turn to Fig. [8] It shows an example similar to Fig. [6J for a different value of Xf, here Xf = 30. The 
phenomenology is similar to that discussed for Fig. [5J the main difference lies in the much lower values of fi^ where 
the isocurvature perturbation is maximal. Since Xf is larger than in Fig.[6j the post-freeze-out dark matter abundance 
is (much) smaller, hence a lesser degree of contamination by curvaton decay is sufficient to produce a substantial 
isocurvature perturbation. One can use the formulae ([3"T]) . (|3"5)l and (|43[) to give a detailed account of Fig. [5J 

Finally, in a third series of figures, we show the same value S^y / Sa\ in the plane of the branching ratio B x and 

initial curvaton abundance • The location of the peak of the isocurvature abundance, as shown in Fig. [9) follows 
the trend described above. Indeed, one sees two main regions where the relative isocurvature perturbation is of order 
unity, one at large branching ratio and small curvaton abundance, which corresponds to the case in which the curvaton 
decay produces a substantial amount of dark matter but negligible radiation, and the other at small branching ratio 
but large curvaton abundance, corresponding to the opposite limit. 

The right panel of Fig.[5]shows a plot similar to that in the left panel, albeit for a much larger value r CT /r xx | f = 10~ 3 , 
meaning that curvaton decay occurs shortly after dark matter freeze-out (at Xd ~ 10 2 ). Consequently, secondary dark 
matter annihilations due to curvaton decay are more effective and restrain the part of parameter space where the 
isocurvature perturbation is maximal. The lower region is also pushed to larger values of Q$ since, at given Cl„ , an 
earlier curvaton decay means a smaller value of f2^ d . 
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FIG. 9: Contour plot of S^/Sjj\ in the plane B x , for xi = 20. Each contour represents a decrement by a factor 2. In the 
left panel, V a /T xx ^ { = 10" 6 ; in the right panel, r CT /F xx | f = 10" 3 . 

V. CONCLUSION 

In this article, we have constructed a comprehensive model of curvaton cosmology including radiation, WIMP 
dark matter (neutralinos) and a massive curvaton field modeled as a perfect fluid with a dust-like equation of state. 
We have computed the evolution of cosmological fluctuations through dark matter freeze-out and curvaton decay by 
approximate analytical methods as well as with a fully exact numerical treatment. As a result, we determine the final 
curvature and isocurvature perturbations and define the parameters that control their amplitude. 

First of all, we have confirmed some results that are intuitively clear. There is obviously no production of isocur- 
vature perturbations if the curvaton contribution is completely negligible or if, on the contrary, radiation and dark 
matter both come entirely from curvaton decay. Another case where the final S xl vanishes is when the curvaton de- 
cays before dark matter freeze-out: because dark matter and radiation then share thermal equilibrium after curvaton 
decay, the isocurvature mode is erased in agreement with the theorem of Weinberg on the approach to adiabaticity 
for interacting fluids fill ]. 

Our main finding is that a substantial radiation-dark matter isocurvature mode can be produced if either the dark 
matter or the radiation (but not both at the same time) mostly originates from curvaton decay. Quantitatively, the 
first situation corresponds to B x £l< d > il x d and (1 — B x ) il< d <C £l^ d (we recall that B x denotes the branching ratio 
of curvaton decay into dark matter, and that the symbol <d is intended to mean "immediately before curvaton decay"). 
Despite the fact that the curvaton is underdominant at its decay, its influence on dark matter is important because 
the dark matter energy density is itself even smaller. The second regime corresponds to B x ~ and £!< d > il< d . In 
this case, at decay, the curvaton perturbations are mainly transferred to radiation. The curvaton perturbations also 
modifies the freeze-out at the perturbative level, thereby contaminating the dark matter fluctuations. Nevertheless, 
a significant radiation-dark matter isocurvature mode (as measured relatively to the curvature perturbation) is also 
produced in this limit. 

Another important effect that we have observed and modeled is that regenerated dark matter at curvaton decay may 
annihilate, leading to a partial (or even complete) erasure of any previously existing S xl isocurvature perturbation, 
even if the curvaton energy density at that time is small, il^ d <C 1. Assuming the curvaton decay to be instantaneous, 
we have provided a new analytical formula describing this phenomenon which reproduces fairly well our exact numerical 
estimates. These late time annihilations further restrict the region in the parameter space of the model (in particular 
the value of the initial curvaton energy density) where significant isocurvature perturbations can be produced. 
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APPENDIX A: PERTURBATIONS TO FIRST ORDER 



In this first appendix, one briefly recalls how the gauge-invariant equations of motion of the (scalar) perturbations 
used in the main text, see Eqs. ([T!?|) -- (fl~5]) . can be obtained. The scalar part of the perturbed metric can be expressed 

as mill 

ds 2 = a 2 {f]) {-(1 + 20)d77 2 + 2d l BAx l Ar 1 + [(1 - 2tp) 8 t] + 2d l d J E] dx l dx j } , (Al) 

and depends on four unknown functions: <f>, B, ip and E. In the previous expression, a(rf) is the Friedmann-Lemaitre- 
Robertson- Walker scale factor and 77 denotes the conformal time. The most general (infinitesimal) change of coor- 
dinates (or "gauge-transformation") that can be constructed with scalar functions (given here by £° and £) is given 
by 

r ] ^fi = ri + t {r 1 ,x k ), x { -f x l = x i + S ij d^(rj, x k ) . (A2) 

Calculating the Lie derivative along the four vector defined above, we find that the four scalar functions used to 
construct the scalar perturbed metric transform according to 

$=<t>+f + -e, s = B-e+e, ^=^--e, e=e+z, (as) 

a a 

where a prime denotes derivative with respect to conformal time. As a consequence, the following combinations, 
known as the Bardeen potentials 



- [a (B - E')]' , * = $ - — (B - E') , (A4) 
a a 



are gauge invariant, that is to say, $ = $ and ^ = ty. In particular, $ can be viewed as the relativistic generalization 
of the Newtonian potential. 

One must also construct gauge-invariant combinations for the scalar quantities appearing in the perturbed stress- 
energy tensor. The rule of transformation of the stress-energy tensor is that of any two-rank tensor and has already 
been given above. In the following, we consider the case where several fluids are present and we denote each species 
by the index "(a)". Then, straightforward manipulations lead to: 

Sp(a) = &P(a) + P(a)£° V (a) = v (a) ~ C *P(a) = S P(ct) + P[ a )^ , ( A5 ) 

where pi a -\, P( a ) and Vt a \ are respectively the background energy density, pressure and velocity of the fluid "(of)". For 
the quantities describing the matter, these expressions play the same role as Eqs. (|A3|) for the metric perturbations. 



17 



From these expressions, one can easily construct gauge-invariant density contrast, pressure and velocity perturbations, 
namely 



Ap (a) = 5p (a) + p{ a) (B - E') , A P(a) = Sp (a) + p{ a) (B - E'), V (a) = v (a) + E' . 



(A6) 



The equations of motion for these quantities can be obtained from the perturbed Einstein equations or from the 
time and space components of the conservation equation, the two methods being equivalent thanks to the Bianchi 
identities. 

In the presence of multiple interacting fluids, which is the case of interest in the present article, one must include 
the various energy transfer rates between the different fluids. One may express these transfers in the form: 



M (a) 



?(a) + Y (a) 



(A7) 



where the vector Q" a ) represents the case of a decay and must be linear in energy density while Yj\ describes the 
case of annihilations and must thus be quadratic in energy density. This suggests the following covariant writing for 
the transfer related to decay 



hex) 



-r 



(A8) 



^a(a) being the decay rate of species a to species (a); the above equations holds for a source, not a loss term. In 
addition, the previous expression is valid if the stress-energy tensor is that of a perfect fluid. As far as annihilation 
terms are concerned, the analog y/\ of the transfer rate Qf a ) can be written as 



m 



(«) 



(a)^ (a) (a).cq^ (a),eq 



U/3 



^XX V )(a) 



1)1 



(a) 



2 2 
P(a) ~ P(a),eq 



(A9) 



In this expression T/% oq is the stress-energy tensor of the fluid "(a)" in thermal equilibrium. In particular, the 
corresponding energy density is given by Eq. ([3|). As before, the last equation has been obtained postulating a perfect 
fluid. Of course the total stress-energy tensor is conserved and, therefore, one has: 



Yl [)>) + ^(a) 



(A10) 



One advantage of the previous covariant formalism is that the transformation properties of the vectors Q M and 
can now be evaluated quite easily. Using Eqs. (|A8|I and (jA9|) . one finds: 



SQL) = r CT(Q) [STfuu + T^5u v ) 



(All) 



5Y,\ 

{at) 



( a XX V )(c) 



1(a) 



61, ^xl , \Ur + IT \\ol , \Ua + 1, \\-L i \OUr — oly \ \1, \ ua — l; -, \o±, \ un 

(a) A (a) P (a) (a) " i a ) ( a ) P (°0> ec l ( a )i cc l ^ ( a )> c q ( a ): e q " 



-iA/3 



At) 



-Tf . x tM 5u 

(a),eq A (a).cq P 



(A12) 



Working out the above expressions for the time component, this gives 

(°xx w )(a) 



8Q(a)0 = a^a(a) [^Pa + P. 



$ Y (a)0 



|-2a [/?( a )5p( a ) - P(Q), cq ^(Q),cq] - a 



2 2 
P(a) _ ^(a),eq 



(A13) 

where we have used that i5m° = —(f>/a. Then, the time component of the conservation equation leads to the following 
gauge-invariant equation, for instance for the dark matter component: 



Ap' + m (Ap x + Ap x )-3 (p x + Px ) + { Px + Px ) O k d k V x = aT x (Ap x + p x $) 



^f^a [2p x A Px - 2p XtCq Ap XtCq + (p 2 x - p Xt J $] 



(A14) 



On large scales the gradient of the velocity term can be neglected. Then, the above equation exactly reduces to 
Eq. (fT2"|) used in the main body of the text. Other equations (fT5)l and (fH|) are obtained along the same lines. 
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APPENDIX B: JUNCTION CONDITIONS 



The junction conditions allow us to derive how the "conserved" quantities cross freeze-out. Assuming that the 
spatial transition hypersurface E is defined by the condition that some function q = q + Sq remains constant, its 
normal n M is given by = d^q and, then, two conditions must be fulfilled [20, 21] so that the two space-time manifolds 
along £ can be joined without a surface layer: the induced spatial metric _Uj obtained from A-p V = + n^n v and 
the extrinsic curvature (or second fundamental form) K^j should be continuous on £, i.e. 



The extrinsic curvature is defined as: 



K„ 



-pv v 



V v n") 



(Bl) 



(B2) 



where C n denotes the Lie derivative with respect to the normal . In order to compute the system of coordinates 
(i.e. the gauge) and the vector (i.e. the surface of transition) have to be specified. Different choices for ra M lead to 
inequivalent junction conditions. The previous equations are the general rules that must applied in order to find the 
matching conditions. Expressing the condition [6 _Ly] , = for diagonal and off-diagonal terms leads to: 



hH»k)]± = 



(?7o)=0, [E(v s )] ± = [E( VQ )} ± =Q, 



(B3) 



where we have used that r/^, = rjo + St). Repeating the same calculation for the second fundamental form gives two 
other (gauge-invariant) junction conditions, namely 



V/ + h<p+(h' -n 2 ) 5 4 



B-E'+ 6 4 
%. 



M = o. 



(B4) 



The equations (|B3|) and (|B4[) represent the complete set of matching conditions for density perturbations. Let us now 
derive the consequences of these equations. Since [ip]± = —H [Sq/q' ] ± , [B — E'] ± = — [Sq/q' ]± and [\l/] ± = [$] ± = 
[ip]± — Ti [B — E'] ± (the first equation of this last relation simply comes from the fact that we assume <£> = "J), we 
deduce that 



[*]- 



0. 



(B5) 



This condition does not depend on which quantity q we choose and is therefore hypersurface independent. On the 
other hand, to obtain the other condition, we must specify this quantity. For the particular case of dark matter 
annihilation freeze-out, the transition is determined by the quantity 



q ° = H- 



Before freeze out one has r/H oc T a /p^ 2 and this leads to 



-ip-H 



Sq 



C 7 <f + 



<f 



30. (C< f -C 7 <f ) 
[2 (2 - a) - n a ] 



o(n, 



After freeze-out, one has Tn x (a 



x\ u xx 



-ip-H— 
% 



r) x p x . This leads to 
2 



C >f + 



1 



2--n c 



^x 



-o CT e f -2(i-n CT )C 7 >f 



2 + n t 

Equating these two expressions and using the fact that £5* = and £< f — £> f , one arrives at 

(«-3)0 CT , f 



/■>{ — /-<t 



2 (a - 2) + Q, 



erlf 



(B6) 



(B7) 



(B8) 



(B9) 



where the background quantities are evaluated at the time of freeze-out and where we recall that a = m/T + 3/2 
evaluated at the freeze-out. We have thus recovered, by means of the junction conditions, Eq. (f2"2")l used in the text 
and derived for the first time (although using different techniques) in Ref. Q. 
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APPENDIX C: APPROACH TO ADIABATICITY FOR INTERACTING FLUIDS 



In this last appendix, we discuss the theorem proved by Weinberg [111 ] , according to which any pre-existing isocur- 
vature perturbation between two interacting fluids is rapidly erased at thermal equilibrium. In addition, we explicitly 
calculate the typical time scale of this phenomenon in the present context and compare the analytical prediction to 
numerical calculations. 

The number density of dark matter particles evolves according to n' + ZTin = Y /m x , with Y /m x — 
—a(<r xx v) (n 2 — rip q ) , in the absence of curvaton decays but accounting for x~ X annihilations. This equation cannot 
be solved explicitly but useful information can be obtained by linearizing it. For this purpose, one writes n = n cq + Sn 
where Sn <C n cq . Then, the equation obeyed by Sn reads, to first order in Sn/n cq : 



dt 



(a 3 (5n) + 2{a xx v)n eq a 3 Sn 



dt 



(a 3 n cq ) , 



whose solution can be expressed as: 

Sn = Sn(ti) 
where we have defined the function: 



a(ti) 



1 3 



a(t) 



n(U: t)-— I d9— (a 3 n cq ) n(9; t) 



Tj(t'; t) = exp (~ 2 J dr ( a xx v ) n <*ij 



(CI) 



(C2) 



(C3) 



As long as x < xr, the integrand (annihilation rate) (a xx v)n eq is very large compared to an inverse Hubble time, 
hence the function n(t'; t) collapses to zero on a very short timescale ~ \{a xx v)n cq ] 1 . One can thus Taylor expand 
the integral to first order and write: 



r)(t'; t) ~ exp [—{t — t') x 2(a xx v)n c 



(a xx v)n eq {f) 



(C4) 



using the following representation of the Dirac function, S t (x) — e~' :E '/ c /(2e). 

The first term on the r.h.s. of Eq. ()C2|) corresponds to the difference at t; between n eq (ti) and the solution of the 
evolution equation for n. This term collapses on the above short timescale due to the presence of the function rj(ti; t). 
This timescale can be written in a compact way in terms of the e— fold number N, using the change of variables 
d< = dN/H: 



t) 



exp 



N-Ni 
AN 



with: 



AN = 



-Xf-\-X\ 



2{a xx v) ti n eq (t i ) 2^/xf2 



(C5) 



(C6) 



For Xi = 3 and Xf = 20, the number of e— folds necessary to erase any initial departure from the solution n is: 
A7V ~ 2.67 x 10 -9 , very small indeed as compared to unity. 

The second term on the r.h.s. of Eq. (|C2[) represents the difference between n and n oq . It will prove interesting for 
the discussion that follows to estimate the relative magnitude of this term. Assuming that 5n; = 0, one can re-express 
this term by changing variables dt = d\n(x)/H in Eq. (|C2[) : 



5n{x) 



dx'l- 



3/2 



n eq (x')n(x' ; x) 



1 x - 3/2 fx 

2 x~ 



1/2 



— I e~ 

Xf 



^ , eq(^) ) 



(C7) 



where we have approximated r/(x'; x) with the delta function in the last equality (the factor of 1/2 comes from the 
fact that the Dirac function is peaked over the upper boundary of the integral). Hence Sn/n eq ~ O (e x ~~ Xf ) <C 1 when 

X < Xf. 

Finally, by integrating by parts Eq. (|C2|) , one can write the global solution n as: 



n{t) 



n(ti) 



r 7 (t i ;t) + -( ( 7 xx «> / a 3 (9)n 2 cq (9)n(9;t)d9 



(C8) 
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Here as well, the first term on the r.h.s. represents the initial value which collapses on a time scale (2(a xx v)n eq )~ 1 , 
while the second term can be shown to approximate n cq to within the above Sn. 

Let us now turn the perturbed quantities. At the perturbed level the equation V^n^ 1 = gives 



d Sn° + diSrf + iHSn n + (-3tp' + d.d'E' + </>') n° = . 



(C9) 



The quadrivector n M can be written as n^ — nu^ and, therefore, one has Suq = —a (Sn + n</f>). Then, one can define 
the following gauge-invariant quantity 



An = Sn + n' (B - E') , 



(CIO) 



and, on large scales, Eq. (|C9|) reads 



An' + 3V. An - 3n*' = . 



(Cll) 



This equation is nothing but Eq. (|Al4jl for a pressureless fluid since then p = m x n. Again if one takes into account 
the interaction term, then the above equation becomes 



SY 

An' + ZHAn - 3n$' = — , 



(C12) 



where (using again the relation p = m x n) the explicit form of the annihilation term is given by Eq. (j A13|) . namely 



= —2a((T xx v) (nAn — n CQ An CQ ) H $ , 



(C13) 



■■eq 



n cq A cq , and A Gq has been defined in Eq. f|16[). In particular, since An eq is proportional to Ap 7 , its 



(C14) 



where A? 

gauge- invariance is automatically guaranteed. In Ref. [9j, the following quantity 

An 

C x ee -$ - n— , 



has been defined. This quantity is obviously gauge-invariant and, if Q M and are not considered, is simply equal to 
C x . Then, let us consider the following combination 



_d_ 

dij 



A^ 

P'-r 



n' + 3Hn - 



SY 

An' + 3HAn - 3n$' = . 



(C15) 



Working out this expression, one obtains 
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drj 
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= -3H An - n 



— 2a(a xx v)nAn + 2a(a xx v)n eq An cq 
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m x p> 7 
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$' + n 



p\ 



n 



Ap 7 y 
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(C16) 



This formula is exact. At this stage, one should make some approximations adapted to the situation at hands. Let 
us first assume that we are before freeze-out. According to the previous considerations, one can neglect Y since its 
magnitude is of order Sn/n eq « 1 (i < Xf) as compared to the others. The previous equation then reduces to 



dn 



^An - n '~T^j = (-3?* - 2a(a xx v)n cq ) ^An - n'^y 1 ^ + 3" 



' A P"f , v ( A P-y \ ' 



P-y \ Pi 



(C17) 



Now, we further assume that freeze-out occurs in a radiation-dominated era. In this case, the quantity defined 
by £7 = — $ — TLAp^l p'^ is conserved hence the term proportional to 3n vanishes. Therefore, using the fact that 
n' ~ —Tin (x + 3/2), the final equation reads 



_d_ 

d'i] 



. An n' A/3 7 
In I f- 

n n p' 



= [ x - - I TC - 2a{a xx v)n eq . 



(C18) 
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N (e-folds) 

FIG. 10: Evolution of the isocurvature perturbation S XJ [relative to the initial S^-y] during the first fractions of e— folds at 
x ~ 3. The thick solid line shows the results of the numerical integration, and the dashed line that of the analytical calculation, 
given by Eq. (|C19|I . 

However, before freeze-out, the first term in the right hand side of the above relation is much smaller than the second 
one and, hence, can be safely neglected. Finally, one arrives at 

4 - Cr = (4 - C 7 ). fqrl^to *) . ( C19 ) 

where the function r)(ti\ t) has already been computed in Eqs. (|C4j) and (|C5|) . 

This analytical result is compared to the numerical calculations in Fig. |TD] Indeed, the agreement is excellent, 
and the numerical calculations confirm the erasure of any initial isocurvature perturbation S xl on the e— fold scale 
AN - 3 x 1(T 9 computed in Eq. JC6J. 



